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A simple model of ID structure based on a Fibonacci sequence with 
variable atomic spacings is proposed. The model allows for observation 
of the continuous transition between periodic and non-periodic diffraction 
patterns. The diffraction patterns are calculated analytically both using 
"cut and project" and "average unit cell" method, taking advantage of the 
physical space properties of the structure. 
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1. Introduction 

For nearly 100 years, the analysis of diffraction patterns of solids has 
been an essential tool for studying solids, since the diffraction pattern of 
a solid is essentially the (squared) Fourier transform of the set of atomic 
positions. Classical crystallography considered periodic structures, whose 
diffraction patterns consist entirely of sharp Bragg peaks. The Fourier trans- 
form of such a periodic set can be computed from the relevant unit cell. The 
discovery of quasicrystals showed that discrete diffraction patterns 
sociated not only with periodic structures but also with a large family of 
solids that have no discrete translation symmetry - quasicrystals. This fact 
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was incorporated into a new definition of "crystal" proposed in 1992 by the 
Commission on Aperiodic Crystals established by the International Union 
of Crystallography: a crystal is defined to be any solid with an essentially 
discrete diffraction pattern. 

Diffraction patterns for periodic and aperiodic crystals differ in substan- 
tial ways - for instance, the diffraction patterns of quasicrystals may exhibit 
"forbidden symmetry" . It is therefore illuminating to consider a model that 
interpolates between periodic and aperiodic structures, and observe how the 
diffraction pattern changes. 

We consider such a one-parameter family of structures in this pa- 
per. Specifically, we consider a fixed (Fibonacci) sequence of two types 
of "atoms" , and vary the amount of space around each type of atom, while 
keeping the overall density fixed. The control parameter k is the ratio of the 
two allowed distances between nearest neighbors. In all cases, the diffraction 
pattern is discrete, and the locations of the Bragg peaks are independent 
of K. However, the intensities of the peaks are K-dependent. When k is 
rational, the intensities form a periodic pattern, while when k is irrational, 
the diffraction pattern is aperiodic. We compute this diffraction pattern in 
two independent but equivalent ways: a) by recovering periodicity going to 
higher dimension (the "cut and project method" deeply discussed in many 
papers: [117113112111311111131131171111112]); b)using the concept of the 
reference lattice. 

These results are in accordance with the ergodic theory of tiling spaces. 
It is known that the Bragg peaks of a tiling T occur at eigenvalues of the 
generator of translations on the hull of T (i.e., the space of all tilings in the 
same local isomorphism class as T) 8 . It is also known 24^ that the hulls of 
modified Fibonacci chains with the same average spacing are topologically 
conjugate, hence that their generators of translations have the same spectral 
decomposition. The question of when and how such a modification affects 
the dynamical spectrum was addressed for one dimensional patterns in [51, 
and for higher dimensional patterns in 6 . (It should be noted that for a 
substitution tiling whose substitution matrix has two of more eigenvalues 
greater than 1, a generic change in tile length will destroy the Bragg peaks 
altogether, in sharp contrast to the behavior of modified Fibonacci chains, 
other Pisot substitutions, and other Sturmian sequences.) 

Ergodic theory says nothing, however, about the intensities of the Bragg 
peaks. Although the spectrum of the generator of translations is compli- 
cated, for special values of the control parameter some of the peaks may have 
intensity zero, resulting in a simpler diffraction pattern. The calculations 
in this paper demonstrate that this does in fact happen. 
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2. The Modified Fibonacci Chain 

The properties of Fibonacci sequences have akeady been thoroughly 
studied (see e.g., [2S1)- They are sequences of two elements A and B ob- 
tained from a substitution rule: 

A — > AB; B — > A. (1) 

Let 

Pm — {Pm^Pm) population vector, where tells how many 

elements of type X are among the first m terms of the sequence. Of course, 
Pm + Pm — There are an uncountably infinite number of Fibonacci 
sequences, but all have the same local properties and the same diffraction 
pattern. It is easy to see that every Fibonacci sequence has 



lim ^ = T, where r = — ^ (2) 



For definiteness, we will work with the sequence (HSl) 



m 



T 



m 



T 



Pm = i — ,m- — ). (3) 



Here ||-|| is the nearest integer function: If m G Z and m < x < m+ 1 then: 

m if xS[m, m + 1/2), 

m+1 if x G [m + l/2,m + 1). 

Now pick two positive numbers (also called A and B) that determine 
the space between each or "i?" atom and its predecessor. That is, the 
atomic positions are given by 

Xm=Pm-{A,B). (5) 

We call the sequence {xm} a modified Fibonacci chain. If A = B, then the 
atoms are equally spaced, and this is simply a periodic array. If A/B = r, 
then the atomic positions are those obtained from the canonical "cut and 
project" method. By varying the ratio A/B, we interpolate between these 
two cases. Note that the average atomic spacing is We will keep this 

average spacing fixed and consider the one-parameter family 

A = T--, B = l + €, (6) 
r 



depending on the control parameter e. 
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3. 2D Analysis of the Modified Fibonacci Chain 

The modified Fibonacci chain can be obtained by a "cut and project" 
method with a nonstandard projection. From this construction we can 
compute the diffraction pattern. 

Let be the 2 dimensional square lattice with spacing u. The Voronoi 
cell of each lattice point is a square. Let /q be the line y = x/t through the 
origin, making an angle /3o = cot~^(T) with the x-axis. Let X be the subset 
of whose Voronoi cells are cut by Iq. It is well known (see e.g., US]) that 
X = {x \ X = lypmjm e Z} where Pm are population vectors given by (jSJ. 

Let la be the line through the origin making an angle 

/? = /3o + a (7) 

with the X-axis, and define to be the orthogonal projection onto la- 
Finally, let A be the projection of X onto la- 

A = n„(x). (8) 

The set A is then a modified Fibonacci chain; the procedure was presented 
in the figure ^ The two distances are 

A = ucosl3, B = u sin /3. (9) 

Their sequence is fully determined by X and does not depend on a. The 
average distance between nearest neighbors in A is 

1- • ^) ^/^^ nn\ 
hm — — = v cos a. (lUj 

m^oo -(1,1) r + 1 

To keep this average spacing constant we take 



(H) 

COS a 

The angle a, the displacement parameter e of @ and the ratio k = A/B 
are related by 

'^"^ Ttan(a), (12) 



K + r - 1 

r — tan(a) r + e(l — r) , , 

K = cot /?) =- V4 = 7^ 13 

^'^^ 1 + Ttan(a) 1 + e ^ ' 

tan(a) =- = . (14) 

T KT + 1 
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We assume that every point of the set A is an atom with scattering 
power equal to unity. Our aim is to calculate the diffraction pattern of such 
a structure. We begin by calculating the 2-dimensional diffraction pattern 
of X. The diffraction pattern of A is then a section of the diffraction pattern 
of X along the direction 1^ (figure |2) . 

To get the diffraction pattern of X we note that X is ^2 times the 
characteristic function of a strip of width 

T + 1 , , 

K = 15 

COS a 

around /q- The Fourier transform of a product is the convolution of the 
Fourier transforms, and the Fourier transform of a lattice is the reciprocal 
lattice. The diffraction pattern of X in 2D has normalized intensity 

Tn\ V- /'sin((/i<^|fc - A?„^^ |)/2)n2 n / 1^^ 

nix my i^^rnxniy])/ ^) 

where km^my = ^{^x,^y) and rux, my G Z label points of the reciprocal 
lattice to Along the direction la, the peaks can be observed at positions 

fjPhys _ ^rng^my ' (t, 1) _ ^^ T^x + my 

cosaV7T2 r + 2 ■ 

Note that these positions are independent of a (or e or k). Their intensities 
are 

/m,m, = ( — ), W= ^ „ r—^ • (18) 



sm.wV'^ /iQ,fcm^m,„ • (— 1 — rtana,r — tana) 

> w = , 

w J 2\/fT2 



After simplifying and rewriting in terms of e we obtain 

r + 2 



4. Structure Factors and Average Unit Cells 

The concept of a reference lattice has previously been proposed in 
Suppose we have a 1 dimensional Delone set {r^}. Its points represent 
atoms, whose scattering instensities are equal to unity. With some appro- 
priate assumptions on the sequence we get the following expression for 
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the structure factor: 



1 ^ 1 ^ 

F{k) = lim — y'exp(iA;r„) = lim — : exp(iA;n„) 



n=l 



n=l 



A/2 
-A/2 



P(n) exp(ifcu)dii, 



(20) 



where P{u) is the probabihty distribution of distances Un from the atoms 
to reference lattice positions mX, A = 2ii/k. That is, 



A. 



(21) 



We call the series u„ the displacements sequence of r„ (induced by the 
reference lattice with period A). Any series u'^ such that 



A 



A, 



(22) 



will be called an unreduced displacements sequence (of r„). 



Theorem 1. Let r„ = a„ + &e a sum of two real series. If dn is a 
displacements sequence of a„ induced by a reference lattice, then 

u'n = dn + f3n (23) 

is an unreduced displacements sequence of rn induced by the same lattice. 

Proof. Let A be the period of the reference lattice. We have to show that 



'•n — '71 



A 



A = n' 



A 



A. 



(24) 



Note that for any real number x,y, \\x — \\y \\\\ = \\x\\ — ||y||, since ||y|| is 
an integer. We can write the right hand side of (|24|) as 



A 



Pn + OLn- ||an/A|| A 



A 



A = r„ 



which is the left hand side. 



A, (25) 



□ 



The quantity P{u) may be viewed as a probability distribution for an 
average unit cell. The structure factor for the scattering vector k is just the 
first Fourier mode of this distribution. 



buczek printed on February 2, 2008 



7 



Unfortunately, for each scattering vector we get, in principle, a different 
average unit cell and a different distribution. However, the structure factor 
for mk, m G Z can be computed from the reference lattice for k; it is the 
m-th Fourier mode of the distribution P{u). Thus, a single average unit cell 
is sufficient to analyze structures whose scattering occurs at multiples of a 
fixed scattering vector k^. This situation includes, but is not limited to, the 
case where the original point pattern was periodic with period 271 /k^. 

For modulated structures (including quasicrystals), there are usually 
two periods, o and 6, which may be incommensurate. Using two reference 
lattices, the first one having periodicity a and the second having periodicity 
6, the structure factor for the sum of two scattering vectors ko = 271 /a and 
qQ = 27r/6 can be expressed by: 



1 ^ 

-^(^0 + go) = lim — exp(i(/co + qo)xn 



ra=l 
N 



lim ^ exp{i{koUn + QoVn)) 

I irv-i /V ' ^ 



n=l 
a/2 pb/2 



a/2 J-b/2 



P{u,v) exp{i{kQU + qQv))dvdu, (26) 



where u and v are the shortest distances of the atomic position from the 
appropriate points of two reference lattices and P{u, v) is the correspond- 
ing probability distribution, which thus describes a two dimensional av- 
erage unit cell. Likewise, the structure factor for a linear combination 
n/co + Tnqo, n, m S Z is given by the (n,m) Fourier mode of P{u,v). This 
means that the average unit cell, calculated for the wave vectors of the main 
structure and its modulation, can be used to calculate the peak intensities 
of any of the main reflections and its satellites of arbitrary order. Using 
and its generalization, it is possible to calculate the intensities of all 
peaks observed in the diffraction patterns. 
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5. ID Analysis of the Modified Fibonacci Chain 

Let a = ^ij, the average spacing between atoms in a modified Fibonacci 
chain, and let b = ra. The positions of atoms in our chain are: 



Pm ■ {A, B) 

m 



A + mB 



T 

{A-B) 



B 



+ mB = {A-B) 



m 



m 1 
- + 2 



m 



+ mB 



{A - B)- +mB + {A - B)( - - \ - + - 

ma + UoM{ma), 



(27) 



where M{x) 



{| + ^} is a periodic function with period b and 



Uo 



r - 1 



er. 



Here [x] = \\x — k\\ is the greatest integer function and {x} 



(28) 
[x] is the 



fractional part of x. 

When Uo ^ 0, the modified Fibonacci chain is thus an incommensurately 
modulated structure. A similar derivation for e = can be found in [2^. 

We are going to construct a two dimensional average unit cell based on 
the two natural periodicities for Xm- o, and h. For the scattering vector 
fco = 27r/a, it is obvious that the series u'^ = UoM{ma) is an unreduced 
displacements sequence of Xm induced by the reference lattice with period 
a. 

Next we consider the (one dimensional) average unit cell for the scatter- 
ing vector go = 27r/6. The series 



ma 



ma 



-bM{ma) = —u'^^ 



(29) 



is an (unreduced) displacement sequence of ma induced by a reference lattice 
with period b. Using Theorem ^ we immediately get that the series 



'^'m = + ""m = (^0 — b)M{ma) = VoM{ma) = ^-u^ 



(30) 



is an unreduced displacements sequence of Xm induced by a lattice with 
period 6, where 



i = Vo/uo 



l-eT2 ■ 



(31) 
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By Kronecker's theorem (see the series u'^ is uniformly distributed 
in the interval [— \uo\ /2, \uo\ /2]. As pointed out by Elser (,9|, for a more 
precise discussion see also |25| ) the uniformity of this distribution is crucial 
for our deliberations. Likewise, the series v'^ is uniformly distributed in the 
interval [— \vq\ /2, |i;o| /2]. 

The structure factor (see (UHl) is 



/a/2 i'b/2 
/ P{u,v)eyiY){i{nikQU + n2qQv))(lv(lu. (32) 
-a/2 J-b/2 



The unreduced displacements sequences u'^ and v'^ can be used to calculate 
P{u,v). However, this cannot be done directly because their terms may lie 
outside the average unit cell (i.e.: |iio| > a or |wo| > 6). Such a situation is 
shown in figure El We have to reduce the series to the interior of the cell. 
The probability function P{u, v) is nonzero only along segments with slope 
^ (as a result of the strong correlation between u'^ and u'^ given by (jSOI)) 
and has constant value. This last fact follows from the uniformity of the 
marginal distributions. 

The formula H32|) is invariant under the changes 

n— >u + 7ia, u— >f + 726, (33) 

where 71^2 sue arbitrary integers. Likewise, the formula does not change if we 
use P'{u, v) = C5{v — (,u) instead of P{u, v) and we change the area of inte- 
gration from [-a/2, a/2], [-6/2, 6/2] to [-11^01/2,1^01/2], [- \vo\ /2,\vo\ /2]. 
That is, we are free to integrate a part of distribution in the neighboring 
unit cells. 

For integers ni,n2 we compute the location i^„j^„2, the structure factor 
F{Kni,n2) s-iid the normalized intensity / of the corresponding peak: 



27r(rni+n2) r^/r^ sin(w) 

Kmn2 = {niko+n2qo) = F(A„ = > ^ = 

ra w 



where 

w = {niko+n2qoOuo/2 = {Kni^n2-n2qi)uo/2, qi = qo{l-(,) 



2 

(34) 



27rr 



1 - er2 ■ 
(35) 

The integers ni and n2 label the main reflection and its satellites, re- 
spectively. Equations 1^5)) can be used to calculate the positions and 
intensities of all peaks. 

The correspondence with the previous 2 dimensional calculation is given 

by 

ni = my, n2 = rrix — ruy. (36) 
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By equations (|17|) and . the peaks are located at 

_ _ TUyT + m,-myT + l _ m,T + my _ f,^, 

J-'^ni,n2 — ''■^my,mx-my — ^'l ^ -j- _|_ 2 ~ T + 2 ~ 

(37) 

Likewise, substituting into and simplifying yields (fT^ . It must 
be noted that the 2D approach presented here is not new and has been 
already used to describe the transformation between quasiperiodic and pe- 
riodic structures. Please refer to [20] and j28j . 



6. Discussion of the Results 

It has been shown that the deformation rule in physical space changes 
only the amplitude of modulation (equation (E?))). Positions of peaks do 
not depend on the parameter e; only their intensities vary. Using equations 
pi]) and we can easily build envelope functions, which go through the 
satellite reflections of the same order (indexed by n2). The shift of the 
envelope functions is qi, as given by 

The set of positions of Bragg peaks is always periodic, since the spec- 
trum of a one-dimensional dynamical system is an Abelian group. By a 
commensurate diffraction pattern we mean a pattern in which the ampli- 
tudes are periodic as well. However, aside from the special case A = B, the 
Bragg peaks are described by two incommensurate periods, and should not 
be confused with the diffraction of a periodic crystal. For our diffraction 
patterns, one period (of length qi) is connected with envelope functions, 
while the second, with period ko, is associated with peaks ascribed to each 
envelope function. This behavior is characteristic of modulated crystals 
and was discussed in [SOj. Only the first periodicity could assure equality 
of intensities. 

It is convenient to describe our results in terms of the ratio k = A/B. 
When K is rational, say equal to p/q, then every atomic location Xm is 
a multiple of A/p = B/q. Plane waves whose frequencies are multiples 
of 2Trp/A have value 1 at each atomic position, and the entire diffraction 
pattern is periodic with period 2ttp/A. 

Conversely, if the diffraction pattern is periodic, then the underlying 
periods ko and qi must be commensurate. A simple algebraic calculation 
then shows that k must be rational. Thus, the pattern is commensurate if 
and only if k is rational, which corresponds to projecting the set X onto a 
rational direction (k = cot/3; see ©). 

Figure shows the diffraction pattern of an unmodified Fibonacci 
chain. The pattern is clearly non-periodic, as A/B equals to r. For e = 
our approach is identical with that presented in [^. Figure shows the 
pattern for e = —0.7 (corresponding to the average unit cell presented in 
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figure inj. As we can see, analytical calculations of envelope functions are in 
full agreement with numerical calculations of the diffraction pattern. 

Figures IS] and ini show diffractions patterns for k equal to 2/3 and 3/2 
respectively. The regular series of peaks are clearly visible, but the diffrac- 
tion patterns are still quasicrystaline. It is significant that for any value of 
e except 1 — l/r (discussed below) the structure is not periodic in physical 
space, but may have periodic diffraction patterns. 

For e = 1 — l/r one gets fully periodic structure with k = 1, hence 
A = B. Our deformed Fibonacci chain is then simply a lattice, and its 
diffraction pattern is the reciprocal lattice, with period /cq. At this special 
value of K, all the other Bragg peaks have intensity zero. 

It must be also noted that the amplitude of each peak is a continuous 
function of k. In fact, it is infinitely differentiable. As k is varied, there is 
no phase transition between commensurate and incommensurate diffraction 
patterns; the evolution is smooth. As such, with measurement apparatus 
of fixed accuracy, it is impossible to determine whether a given pattern is 
precisely periodic. 

Finally is may be noted that our analysis of the transition quasiperiodic- 
periodic (commensurate-incommensurate) is based on the explicit knowl- 
edge of the structure factor. It has been already shown (JHI) that the 
method advertised (average unit cell) can give the factor for higher dimen- 
sional structures and therefore there are (in general) no obstacles to repeat 
similar analysis in the latter case. Such calculations have not been under- 
taken so far but it seems that they might be based on analytical expres- 
sions for the coordinates of quasiperiodic lattices, derived from periodic or 
quasiperiodic grids (as given for example in |2S|)- 
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Fig. 2. Diffraction pattern of the modified Fibonacci structure is a section of the 
diffraction pattern of X through direction 1^. hP'^y^ = OP is the position of the 
peak, KP determines the intensity. 
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Fig. 3. (a) shows two parameter average unit cell. The distribution P{u,v) is 
non zero only along the thick lines and has constant value. Projections of it onto 
directions u and v determine the probability distributions for scattering vectors ko 
and qo ((b) and (c) respectively), (d) presents set of vectors {{u'jn,v!^) | m G Z}. 
It may happen its elements lie outside the average unit cell and have to be reduced 
to its interior (like for the presented example with e = —0.7). Invariance under the 
substitution H33|l assures that integration of functions on (a) and (d) gives the same 
results. Our parameter space has discrete translational symmetry like a periodic 
crystal. 
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Fig. 4. Examples of diffraction patterns: (a) unmodified Fibonacci chain (k = t and 
e = 0); (b) K ~ 6.836 (e = —0.7). Broken lines present envelope functions. Given 
envelope function goes through satellite peaks of the same order. All envelopes 
have the same shape; their shift is gi. As we can see analytical results are in full 
compatibility with numerical calculations. 
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Fig. 5. Modified Fibonacci chain for k = 2/3 (e w 0.741, qi w 10.816) 
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Fig. 6. Modified Fibonacci chain for k = 3/2 (e w 0.056, qi w 11.913) 



